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Abstract 

The two-particle irreducible (2PI) effective action theories are employed to 
study the strongly fluctuating electron systems, under the formalism of the 
two-dimensional Hubbard model. We obtain the corresponding quantum 
2PI effective action after the original classic action of the Hubbard model is 
bosonized. In our actual calculations, the 2PI effective action is expanded 
to three loops, in which the leading order (LO) and next-to-leading order 
(NLO) quantum fluctuations are included. Numerical calculations indicate 
that the NLO fluctuations should not be neglected when the Coulomb on- 
site repulsion energy is larger than two times the nearest-neighbor hopping 
energy. 
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1. Introduction 

Since the high-temperature superconductivity was discovered in 1986 [if, 
It has been believed that the appropriate model to describe the strongly 
correlated electron systems is the nearly half-filled two-dimensional Hubbard 
model with moderately large repulsion energy U and antiferromagnetic ex- 
change constant J = At 2 fU where t is the site hopping [2j. The Hubbard 
model is just composed of two terms: one is the site-hopping term which 
forms the band structure, and the other is the Coulomb repulsion term that 
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represents the interaction between electrons. Therefore, from its appearance, 
it looks like that the two-dimensional Hubbard model is an simple model and 
it is easy to solved. However, the actual situations are quite different. The 
Hubbard model describes a many-body electron system in which the inter- 
acting potential energy and the kinetic energy are comparable. So we can 
not employ the perturbation theory and treat the potential energy or the 
kinetic energy as a perturbation. 

Lots of efforts have been made to solve the two-dimensional Hubbard 
model and many methods have been developed. For example, the quan- 
tum Monte Carlo simulations the self-consistent approach of conserving 
approximations Q, the variational cluster perturbation theory (if, the func- 
tional renormalization group approach jfjj], and so on. In recent years, an- 
other nonperturbative approach, known as the two-particle irreducible (2PI) 
effective action theory first introduced in the field theory [7|, has attracted 
lots of attentions. The 2PI effective action theory resums certain classes of 
diagrams to infinite order, so nonperturbative effects are included in this 
approach. Furthermore, in the 2PI formalism, the effective action can be 
expanded according to the order of the loop or 1/N in the O(N) model. 
Therefore, it is easy to investigate the effects of the high order contributions 
in the 2PI effective action theory. In the studies of field theories, it has been 
found that the 2PI effective action theory is very successful in describing 
equilibrium thermodynamics, and also the quantum dynamics of far from 
equilibrium of quantum fields. The entropy of the quark-gluon plasma ob- 
tained in the 2PI formalism shows very good agreement with lattice data 
for temperatures above twice the transition temperature The poor con- 
vergence problem usually encountered in high-temperature resummed per- 
turbation theory with bosonic fields is also solved in the 2PI effective action 
theory Furthermore, it has been shown that non-equilibrium dynamics 
with subsequent late-time thermalization can be well described in the 2PI 



formalism (see 10J and references therein). The 2PI effective action has also 



been combined with the exact renormalization group to provide efficient non- 
perturbative approximation schemes The shear viscosity in the O(N) 



model has been computed using the 2PI formalism [12j. Specially, we would 



like to emphasize that due to many people's contributions [13|, [14|, |15|, [16 
it has been clear that the 2PI effective action theory can be renormalized, 
which is quite non-trivial for a non-perturbative approach. 

In this work, we will employ the 2PI effective action theory to investigate 
the strong fluctuations of electron systems. We will expand the effective 
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action to three loops and compute the leading order (LO) and the next-to- 
leading order (NLO) contributions to the fermion and boson self-energies. 
Then we will investigate when the importance of the high order quantum 
fluctuations becomes significant with the increase of the Coulomb repulsion 
energy U. The paper is organized as follows. In section [2] we apply the 2PI 
formalism into the Hubbard model and obtain its effective action. In section 
Owe obtain the LO and NLO contributions to the fermion and boson self- 
energies. Numerical results are presented in Sec. HI In section [5] we give our 
summary and conclusions. 

2. 2PI Effective Action Theory 

We begin with the simplest two-dimensional one-band Hubbard model 
which reads 

H = -tJ2 c\ a c j(J + 4^4^' (!) 

(ij)(T i 

where we only consider the nearest-neighbor hopping t and U is the Hubbard 
on-site Coulomb repulsion energy. We rewrite the interaction term as 

4^4^ = -^(c\a z Ci) 2 + ^(ct f c n + cj ; c 4 ) (2) 

for the convenience of calculations below, where a z is the z component of 
the Pauli matrices. The second term on the right hand side of Eq. (TS]) can 
be absorbed in the chemical potential term. Then we arrive at 

H-fiN = -tj2 4^ - 2- I>I<^) 2 - ^ £ 4^- (3) 

(ij)a i ia 

The classic action corresponding to Hamiltonian given above is 

ia {ij)a ia i 

where creating and annihilating operators in Eq. ([3]) are replaced by their 
Grassmann fields. Including quantum and thermal fluctuations, one obtains 
the generating functional, also known as the partition function, which reads 

Z[r]*, V, J) = [ [dc*] [dc] exp ( - / dr \ ^ c* a d T c ia - t ^ c ia c ja 

J ia {ij)a 
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u 



ia i icr 

+ X)j i (-l)^Vc i ] j. (5) 

i ' 

We will employ the Matsubara imaginary-time formalism throughout this 
work and here = 1/T is the inverse of the temperature. An external source 
term for the composite operator (— l)*c|c*Cj is included in Eq. (jSJ). 

Before we continue the calculations, it would be more convenient if the 
bosonization of the action in Eq. (jSJ) is made first. Up to a constant, we have 

exp { - £ dr [ - U - 5>VcO a + M-lY^Ci] } 

-UB i {-l) i 4a z c i + J^jj, (6) 

where a boson field -B is introduced through its functional integral. Sub- 
stituting Eq. (jHJ) into Eq. (jSJ) and neglecting the irrelevant prefactor on the 
right hand side of Eq. (EJ) , one finds 



Z[r}*, i], J] = I [dc*] [dc] [dB] exp [ - I dr [j] c\ a d T c i(J - t 
. £ + £ (|^ 2 - UB^-iycla* 



(7) 



In order to obtain the 2PI effective action for the Hubbard model, we 
need add two-point sources into Eq. (J7|). Then we obtain 

Z[ V *, V , J, M, K] = j [dc*] [dc] [dB] exp { - [J (c*, c) + I (B) + J int (c*, c, B) 

+r]*c + c*r] + JB + -BMB + c*Kc] } , (8) 
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where the two-point external sources are 



l -BMB = \ \ dndT^Biir^M^T^BjiTj), (9) 

ij 

rP 

C*Kc = / dTjdTj } j } j C* a (Tj)K iatj p(Ti, Tj)Cjp(Tj). (10) 

Here a and j3 are spin indices. In Eq. ([8]) we also used the following abbre- 
viated notations: 

I (c*,c) = / dT[22c* a d T c ia -t22c* a c ja - ^^C* Q Q a J, (11) 
Io(B) = [dr^Bl (12) 

rP 

I iat (c*,c,B) = / dry2(-U)B l (-lYcla z c u (13) 
ri*c + c*r] = / dr^2(ri* a c ia + c* ia r] ia ), (14) 

JB = drJ^JiBi. (15) 

Introducing the generating functional for the connected Green functions 

W[ V *, V , J, M,K} = - In Z[ V *, 7], J, M, K], (16) 

it then follows that 

5W 5W 5W 

_ = b c — — = c c — — = -c* c (17) 

5W 1 SW 

SM^ = 2 (5 ^ C + ^ ) ' jK^- = - {c ^ + S ^^ (18) 

where B c , c c , and c* c are the expected values of fields B, c, and c*, respec- 
tively. G and S" are the propagators for boson and fermion fields. 

The 2PI effective action can be obtained from W through the Legendre 
transformation as follows 

5W 5W 5W 

T[c c ,c* c ,B c ,G,S] = W[r]*,i], J,M,K] - Jj — — f|* tt — r] ia - — 

oJi dr)* a dr] ia 



5W T . 5W 
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= w[ti; v ,j,m,k\- JiBt - vl^ta - cZVia 

-\ti[M{B c B c + G)\ + Tr[K(c c c* c + S)}, (19) 

where summations and integrals are assumed for the repeated indices. The 
trace operates in the coordinate and inner spaces. It can be easily proved 
that 



ST 

sr 



-Ji 



MijBj, 



sr 

Sc* c ~ 
5T 



SG, 



--M, 
2 3 



sr 



ss, 



K 



(20) 
(21) 



Equations (|2"0"|) and f l2~Tj) form a set of self-consistent equations which deter- 
mine the field expected values B c , c c , c* c and the propagators G and S, if the 
effective action can be expressed as a functional of these field expected values 
and propagators. Usually, the expected values of fermion field c c and c* c are 
vanishing when their external sources rj* and rj are zero. We will assume 
c c = c* c = in the following calculations and use B in place of B c without 
confusions. It can be shown that the 2PI effective action can be expressed 
as 



T(B,G,S) 



I(B)+ 1 -TT\nG- 1 + ^Tr(G Q 1 G) 



-TrlnS 



-i 



with 



r int (B,G, S) = -Tr[(-£, 
Here we have 1(B) = Iq(B) and 

5 2 I (B) „_i 5 2 Uc\c 



Gn 1 



r 2 (G,s> 



SB 2 



<J-1 



Tr(So 1 S)+r iQt (B,G,S), 

)S] + r 2 (G, s). 

5 2 L mt (c*,c,B) 



5c*5c 



5 c* 5 c 



(22) 
(23) 

(24) 



J[dc*][dc][dB] exp< 


_ 


\BG' x BMfS- l c+h nt (c\c, 


B) 


} 


J[dc*][dc][dB]exp< 




lBG~ 1 B + c*S- 1 c 




\ 



(25) 

2PI 



where E mean is the mean field contribution to the fermion self-energy; T 2 sums 
all 2PI diagrams. The prominent difference between these diagrams and the 
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perturbative ones is that the propagators constituting these diagrams are 
self-consistent ones G and S, not G and S . But the vertices are bare, 
which are same as those in perturbative diagrams. 



3. Self-Energy of Fermion and Boson Fields 



T 2 in Eq. (j25p receives contributions from infinite 2PI diagrams, some of 
which are shown in Fig. [TJ In actual calculations, it is impossible to sum 
all these diagrams. We have to truncate the effective action to certain order 
(But it is possible to resum some kinds of diagrams to infinite order through 
other methods, see Ref. 17J for more details). In this work we will truncate 
the effective action to three loops, i.e., only the first two diagram in Fig. Q] 
are employed in our calculations. 




Figure 1: Two-particle irreducible vacuum diagrams contributing to the effective action, 
where the solid and wavy lines represent fermion and boson propagators, respectively. 

In the following we will call the two-loop and three-loop diagrams as the 
LO and NLO contributions to I^, respectively. Their expressions are 



-iNLO 



1 

/ dr 1 dr 2 ^) H+l2 ^U 2 G Tl i ltT2 i 2 ti(a z S T2 i 2:Tlil a z S Tl i ltT2 i 2 ^6) 

Jo .•,■„•_• 



I .+«2 +13+14 _T 7-4 f~1 f-1 



xtr(<T S T ^i 4<Tl i 1 (T S Tl i ljT2 i 2 (T S T2 i 2 ^ T3 i 3 u 5*7-2 j 3)T4 j 4 ) , (27) 

where the trace tr only operates in spin space. Up to now, we have expressed 
the effective action as a functional of the self-consistent propagators G and 
S. Then one can employ Eq. f|2ip to obtain the self-consistent equations, 
given by 



s- 1 = S^-Z, 
g = Gn — n, 



(28) 
(29) 



7 



where the fermion and boson self-energies are 

5T int _ 8T 2 (G,S) 

2-* r n "mean ~r r n 

ob do 
« S mean + S LO (G,5) + S NLO (G',5), (30) 
n - ,^ int _ JT 2 (G,S) 
8G ~ 8G 
« n LO (G,5) + n NLO (G',5). (31) 

They are depicted in Fig. [2] and Fig. |3J respectively, whose expressions read 

^Tiiioi,T2i2a2 ~~ 5 ] (~^) U Criii,T2i2 Cr ai/3i'5'Tiii/3i,r2i2/32 Cr /32a2' (^2) 

/81/82 

1314 « 3 a 4 fix—Pi 

x G T4 i 4>T2 i 2 a ai ^ i S T1 i 1 [s liT4 i 4a4 <j a4 i 34 S T4 i 4 [s 4jT3 i ;ia3 o' a3 ^ 3 

X ^13/93,7212/32 ^202' (^) 

n^ii.Tijia = — f/ 2 ( — l)* 1+42 tr(cr 2: 5' T1 j liT2 j 2 (T 2 5' T2 j 2iT1 j 1 ), (34) 



n*i,T2J2 

r-/9 



nSJU = / ^3rfr 4 ^(-f/ 4 )G r3J 3 iT4M (-i; 



Xtr((J S ^1^(7 ^7-414, 7212^" ^T2l2,T3i 3 ^ ^73 43,7"! il ) • (35) 




LO = « *> « f « yNLO = 




Figure 2: LO and NLO contributions to the fermion self-energy. 

Up to now, we have worked in coordinate space. In fact, it is more conve- 
nient to employ the momentum lattices, in which numerical calculations are 
easier. The coordinate lattices and the momentum ones are related through 
the Fourier transformation. For example, for the boson and fermion propa- 
gators we have 



G Tiu,T 2 i2 = P 2 , j=r 52 52 G(uikx, uj 2 h) exp(-za;iT 1 + iuj 2 r 2 
V Nr *—r , 

\ ' Wl«l OJ2K2 
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Figure 3: LO and NLO contributions to the boson self-energy. 



+ik 1 R h - ik 2 R i2 ) } (36) 

S Tl i 1 a 1 ,T 2 i2a 2 = ft 2 Z /HT7x 2 '° 2 ( UJl ^ 1 ' W 2^2) 

V ' CJlfel 0)2&2 



x exp(— zcuiTi + 2w 2 t 2 + ik\Ri 1 - ik 2 R i2 ), (37) 

where N is the number of lattice sites, momentum sums are restricted to the 
Brillouin zone, and the Matsubara frequencies are 

j (2rii + 1)ttT (Fermion) , . 

Ui ~ \ 2n t nT (Boson) ' W 



Employing Eqs. fl36|) and fl37|) . we can reexpress the self-energies in Eqs. ([32 
351) in momentum lattices as follow 



V v / oJifei oj 2 A:2 aja^ 



x5„;^(wi - ui, ki - fci - Q; u 2 - uj 2 , k 2 - k 2 - Q) 
x^, Q2 e-^-^ 0+ , (39) 

Y^o M ^ 2 k 2) = /r 6 -L- E EEEE^ 4 

xG(wifci, a^A^G^sfcs, cu 2 fc 2 ) 

x^a'^^^^i - Wl, fci - fci - Q; w 3 , fc 3 )(T^ Q3 

x 5*03,04 (^3 — ^3? k' 3 — k^ — Q; uj' 4 — W4, k' 4 — k4 — Q) 

X(7 a 4 a' 4 S a' 4 ,a' 2 (^4, ^2 ~ U>2, k 2 - k 2 - Q)<7^ Qa , (40) 

n LO (a; 1 ^ 1) cD 2 fc 2 ) = -r 2 — L-^^^tr^^A;!,^)^ 

xS(oj 2 -u) 2 , k 2 -k 2 + Q; Wi-tUi, fci - fci + Q)o" z ,(41) 
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n-°K*^A) = -r«-L- y, EE" 



(VNY 

xG(uj 3 - u' 3 , k 3 - k' 3 - Q; J 4 - u 4 , k' A - k 4 - Q) 
xtr S(uik 1 ,u 4 k 4 )a z S(u' i , k' 4 ;u 2 + tu 2 ,k 2 + k 2 -Q) 
xa z S(co 2 k 2 ,uj 3 k 3 )a z 

xS(u' 3 ,k' 3 ;ui -ui,h - ki + Q)a z , (42) 

where Q = (tt, tt) t (in units of the inverse lattice spacing) is the nesting 
vector in two dimensions. 

In order to simplify our calculations, we assume that the anti-ferromagnetic 
order parameter B in Eq. (I22p is constant. Then the inverse of the fermion 
propagator can be expressed as the following formalism: 

S~l a2 {^iki,uj2h) = 0C(u u ki)8 UlW2 8 hlk2 8 aiol2 +pD(u x , fci)5 wia;2 4 1)fe2+ Q(7^ ia2 . 

(43) 

It then follows that 

5 , aia 2 ( a; l^l) ^2^2) — 0C{UJi, ki)8 u j lUl2 8k 1 k 2 8a 1 a^r(3D(ui, k 1 )8 ullUJ2 8k lj k 2 +Q 'a 1 a2^ 

(44) 

with 

, \ Cjuukj + Q) 

1 u l) C(cj 1 ,k 1 )C(cj h k 1 + Q)-D(cj 1 ,k 1 )D(co 1 ,k l + QY 1 } 
r)( , x -D(a;i,fci) 

Similarly, for the boson field one has 

G~ 1 (u} 1 k 1 ,u}2k 2 ) = PA^k^S^J^, (47) 
G(uiki,U2k 2 ) = l3A~ l (ui,ki)5 UlU2 5 klk2 . (48) 
Substituting Eqs. QHJ and (jUJ) into Eqs. (jg5)) — (H21) . we obtain 

+SE (wi, fcO^jSa^Lfca+oa^oj,, (49) 

^qi L q 2 (^' 1 ^ 1 ' ^2^2) = £+ L (0)1,^1)^0)2 ^jfca^aioa 

+E] NLO (u; 1 , fci)^A,* 2+ Q< a2 , (50) 
n LO (o)ifci,a; 2 fc 2 ) = n 10 ^,^)^*,, (51) 
n NLO (a)i^i,a; 2 fe) = n NLO (a>i, fcO^fflAfc, (52) 
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u 

x 

u 2 



(7(u>i - wi, h-h- Q)e~ l(a ' 1 - UJl)0+ , 

-i(a>i-u;i)0 + 



/iV) 2 



xD(wi — wi, k± — ki — Q)e 



(53) 



(54) 



u)i — &i — &4 — Q)C{u)\ — u>i — u 2 , k\ — ki — k 2 ) 

, - . v Ji - u 2 , ki — k 2 — Q) + C(u>i - cut, ki — k± — Q) 
xZ)(u>i — u>i — ou 2 , k\ — ki — k 2 )D(ui 1 — ou 2 , k\ — k 2 ) 
+D(cu 1 - ui, k± — k± — Q) 

xC(ui — Ui — u> 2 , k\ — ki — k 2 — Q)D(pi — u>2, k\ — k 2 ) 

+D(u 1 - wi, k-i — k-i — Q) 

xD(ui — u)i — uj 2 , k\ — k\ — k 2 — Q) 

xC{u} 1 -U2,ki-k 2 -Q) , (55) 



4/3-1 



U 4 f3 



x 



(\Ztv) UJ1 i tl W2 & 2 





- u 2 , ki - k 2 


- Q) + C(ui - ui, fci - ki - Q) 


xD(Qi 


— ui— uo 2) ki 


-h 


- k 2 )C(cj 1 - u 2 , ki - k 2 ) 




— ui, ki — ki 


-Q) 




xC(u> 1 


u 2 , ki 


-h- 


-k 2 — Q)C{pi - u 2 , ki - k 2 ) 




— u>i, ki — k\ 


- Q) 




xD(Qi 


— ui— u 2 , ki 


-h 


-h-Q) 


xD(Qi 


- tu 2 , ki - k 2 


-Q) 


(56) 
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and 

u h °(Ox,kx) = - 2u2 —, [cKWi-wi.fci-^ + Q) 



£i : |/.'j 

+Z?(u;i, /ci)/}^ — a)i, fcx — fcj) , (57) 

n NLO (w 1 ,fc 1 ) = ^f/ 4 /?- 1 — [^(^-wi + w^fca-fci + fci) 

xC(ui, kx)C{u 2 + ^1, A; 2 + fci - Q)C(u 2 , k 2 ) 

xC(coi — Ox, ki — k\ + Q) 

+A~ 1 (uo 2 + Ox, k 2 -ki + ki)C(u!, hi) 

xC(u 2 + O x , k 2 + k 1 - Q)D(u 2 , k 2 )D{uj x - Ox, k x - k x ) 

+A~ 1 (lu 2 — oji + Ox, k 2 -ki + ki)D(oji, k±) 

xD(uj 2 + Ox, k 2 + ki)C(u 2 , k 2 )C(ui - Ox, k i -k 1 + Q) 

+A~ 1 (uo 2 —ui + u>i, k 2 -k 1 + ki)D(ui, h) 

xD(u 2 + u>i, k 2 + ki)D(u 2 , k 2 )D(uoi - Ox, kx - kx) 

+A~ 1 (u3 2 — ojx + Ox, k 2 — kx + kx — Q)C(ux, h) 

xD(uj 2 + Ox, k 2 + kx)C(u 2 , k 2 )D(uJi - O x , k x - k x ) 

+A~ 1 (lu 2 — ojx + Ox, k 2 — kx + h - Q)C(ux, h) 

xD{u 2 + Ox, k 2 + kx)D(uj 2 , k 2 )C(ujx - Ox, k x — kx + Q) 

+A~ 1 (uo 2 — ux + Ox, k 2 — kx + kx — Q)D(u x , kx) 

xC(uo 2 + Ox, k 2 + ki- Q)C(uj 2 , k 2 )D(ux - Ox, h - kx) 

+A~ 1 (uo 2 — ojx + Ox, k 2 — kx + kx — Q)D(u x , kx) 

xC(co 2 + Ox, k 2 + ki- Q) 

xD(oj 2 , k 2 )C{u x - Ox, kx - h + Q) . (58) 
Furthermore, from Eq. (|24p it is easy to obtain 

GQ l (Uxkx,UJ 2 k 2 ) = PUS^uzShkai ( 59 ) 

(uixkx,u 2 k 2 ) = 0(-iux + £ki)8 0102 ) 

(60) 

(u)xkx,U) 2 k 2 ) = PUBS^uJk-LM+QVl ia a > ( 61 ) 



with 



= — 2t(cos k x + cos k y ) — [A, (62) 
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where the momenta k x and k y are in units of the inverse lattice spacing. Sub- 
stituting Eqs. ( I4"9"j) — and Eqs. ( 155|) — (IBTj) into the self-consistent equa- 
tions Q2BD and (|2SD, we finally get 

A(w,Jfe) = C/ — /3~ 1 Il LO (u, k) - / 9 -1 n NLO (a;, A;), (63) 
C(u,k) = -iuj + i k - P^i: L + (u,k)- ^- 1 l^ LO (u,k), (64) 
D(u,k) = -UB-^Tt ^^)-^^ ^^). (65) 



4. Numerical Results 

Besides equations (1631 — (|65|) . we still need another two equations to per- 
form the numerical calculations. One is the self-consistent equation for the 
antiferromagnetic order parameter B, which can be obtained from the first 
equation in Eqs. (120)) . as given by 

B = -2/T 1 — D{co,k); (66) 

the other one is the equation for the electron density (number per lattice 
site): 

n = -2(3- 1 = y"C(u,k)e iul0+ . (67) 

We employ a N s xN s xN T lattice, where N s corresponds to the lattice number 
in one spatial dimension and we have two-dimensional lattice number N = 
N s x N s ; N T is the temporal lattice number. The calculations are performed 
in the momentum and frequency (k x ,k y ,iu>) space, which are equivalent to 
those in the coordinate space with spacial periodic boundary conditions and 
temporal anti-periodic (periodic) boundary conditions for fermion (boson) 
fields. 

In order to save computing time, we also employ the symmetries of the 
lattice to simplify our calculations. Propagators and self-energies are invari- 
ant under the following symmetry operations: 

k x ^ ky, k x y k x , ky y ky. (68) 

Therefore, we only need to compute one-eighth of lattice points in the Bril- 
louin zone. Furthermore, we also have the time-reversal symmetry as follows 

A{-u,k) = A(u,k)*, (69) 
C(-u,k) = C{u,k)*, (70) 
D(-u,k) = D(u,k)*. (71) 
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For this reason, it is only necessary to calculate the propagators when uj > 0. 

In this work, we want to investigate the importance of higher-order fluc- 
tuations in the 2PI effective action theory. In another word, we would like 
to study whether the expansions of the effective action in order of loops are 
convergent, when the interaction strength U becomes large. In order to ob- 
tain this goal, we perform three different calculations: the first one is just 
the mean field calculation, i.e., neglecting the influence of the fluctuations 
and assuming I^ = 0. The second one is that only the two-loop contribution 
(the first diagram in Fig. [TJ) to T 2 are included. In this calculation, equa- 
tions (163]) — (157)1 constitute a closed system and are solved self-consistently 
through iterations, but with n NLO = E^ LO = S NLO = 0. As for the third 
calculation, it should have been solving this set of equations self-consistently 
with n NLO = S^ LO = S^ LO ^ 0, but such self-consistent calculation need 
lots of computing time and is beyond our computing abilities. Therefore, in 
the third calculation we substitute the results of the self-consistent equations 
which only include the LO contribution to the fluctuations, i.e., the results 
obtained from the second calculation, into equations (j55|) . ([55]) . and (j58|) to 
obtain the NLO self energies. We should emphasize that this calculation is 
not self-consistent, but it is still reasonable to get some information on the 
convergence of the theory from this calculation. 

Equations (!63|) — ( l67j) are iterated successively to search for self-consistent 
solutions. In order to avoid oscillations during the calculations, we also 
employ feedback in the process of iterations, i.e., the updated value is chosen 
to be a weighted average of the new calculated one and the old one in every 
iteration. We use the following criteria to judge whether the convergence of 
iterations is achieved: 
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where the update of the chemical potential /x is realized through a New- 
ton's procedure at a fixed electron density n. Furthermore, we would like to 
emphasize that there is an exponential factor appearing in Eqs. fl53|) . fl54|) . 
and ( I67p . This exponential factor contains an infinitesimal positive constant 
+ . In the numerical calculations, this constant can not be chosen to be very 
small, because that would affect the convergence of the frequency summa- 
tion, since the infinite frequency summation is cut to be finite in numerical 
calculations. Therefore, in our work we choose a relatively small constant, 
and we also check that the numerical results are insensitive to the choice of 
this constant. 




Figure 4: (color online). Antiferromagnetic order parameter B as a function of the 
Coulomb repulsion U in unit of the hopping t. Here the temperature is chosen to be 
T = 0.05£ and the electron density n = 1. The black solid line with squares corresponds 
to the mean field calculations and the red dashed line with circles to the self-consistent 
calculations including LO fluctuations. For the self-consistent calculations, the Brillouin 
zone is discretized into a 10 x 10 lattice and the temporal lattice number is chosen to be 
N T = 64. 



The antiferromagnetic order parameter B as a function of U, obtained in 
the mean field approximations and the self-consistent calculations including 
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Figure 5: (color online). Dependence of the LO boson self-energy LI (w, k x = 0, k y = 0) 
and the inverse of the boson propagator A(lu, k x = 0,k y = 0) on the Matsubara frequencies, 
obtained in the self-consistent calculations with LO fluctuations included. The Coulomb 
repulsion is U — 2t. These results are based on a lattice of 10 x 10 x 64. 



LO fluctuations, are shown in Fig. HJ In the mean field approximations, the 
frequency summation in Eqs. ( 166]) and f )67|) can be performed analytically. 
One can see that B develops a nonvanishing value with the increase of the 
Coulomb repulsion in the mean field approximations, which means that a 
phase transition occurs and a long-range antiferromagnetic state is formed. 
However, in the self-consistent calculations with LO fluctuations included, we 
do not find any phase transition at T = 0.05t, n = 1, and U = 0.5 ~ 2.5t. In 
order to find the reason why there is no phase transition in the self-consistent 
calculations, we show the LO boson self-energy and the inverse of the boson 
propagator A in Fig. We can observe that II LO is positive at low frequen- 
cies, which results in that A is decreased at low frequencies compared to the 
bare Coulomb repulsion U. Therefore, the effective interaction is reduced by 
the fluctuations at low frequencies, which may account for the reason why 
the phase transition becomes difficult in the self-consistent calculations with 
fluctuations included. 

Figure [6] shows the LO and NLO self-energies as functions of the on-site 
repulsion U. The boson self-energies II are real and are presented in the two 
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Figure 6: (color online). Comparison between the LO self-energies and the NLO ones. 
They are depicted as functions of the Coulomb repulsion U. Temperature T — 0.05i and 
electron density n — 1 are chosen here. k x and k y are in unit of the inverse of the lattice 
spacing. These results are based on a lattice of 10 x 10 x 64. 

top panels. In the symmetrical phase where B = 0, £_ defined in Eqs. ( I49p 
and fl50|) is vanishing. Therefore, we only show S + here, whose real and 
imaginary parts are depicted in the two medium panels and the two bottom 
ones of Fig. |6l respectively. From this figure, one can observe that the NLO 
contributions to II and the real part of £+ are less than the LO ones, but as 
the U is larger than about 2t, the NLO contributions can not be neglected as 
the top-right and medium-left panels show. However, for the imaginary part 
of the S + , we can see that the NLO contributions are completely comparable 
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with the LO ones, even when the U is less than 2t. 




co (t) ro (t) to (t) 



Figure 7: (color online). LO and NLO self-energies as functions of the frequencies. Three 
columns correspond to three different values of k x . In the same way, we choose the 
temperature T = 0.05i, electron density n — 1, and the Coulomb repulsion U = 2.5t. The 
lattice is 10 x 10 x 64. 

In order to make the comparison much clearer, we show the spectrum of 
the self-energies in Fig. Here we choose U = 2.5t. For the boson self- 
energies, we find that the NLO contributions are significant and should not 
be neglected, especially at low frequencies as the top-middle and top-right 
panels show. In the same way, the same conclusions are also appropriate 
for the real part of the E + . Furthermore, we also find that the real part of 
the NLO S + approaches zero at high frequencies. For the imaginary part 
of the E + , we confirm that the NLO contributions are comparable with the 
LO ones, which is clear at the whole region of the frequencies as the bottom 
three panels show. 

Figure [S] depicts the self-energies at several values of the temperature. 
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Figure 8: (color online). LO and NLO self-energies as functions of the frequencies at three 
different values of the temperature. We choose k x = 3ir/5, k y = 0, U = It, and n = 1 in 
all these calculations. These results are based on a lattice of 10 x 10 x 64. 



Here we choose the Coulomb repulsion U = t. We have demonstrated above 
that at this value of U, the NLO contributions to the boson self-energy II 
and the real part of the fermion self-energy S + are quite smaller than the 
LO ones, which are confirmed at other values of temperature as the first two 
rows of Fig. [8] show. In the same way, at other values of temperature, we 
also find that the NLO contributions to the imaginary part of the fermion 
self-energy are comparable to those of LO. 

5. Summary and Conclusions 

In this work, we have employed the 2PI effective action theories to study 
the strongly fluctuating electron systems, under the formalism of the two- 
dimensional Hubbard model. We first bosonize the original classic action of 
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the Hubbard model, then obtain the corresponding quantum 2PI effective 
action. In our actual calculations, the 2PI effective action is expanded to 
three loops. Therefore, LO and NLO quantum fluctuations are included in 
our approaches. We also obtain the LO and NLO contributions to the fermion 
and boson self-energies. They are expressed in the momentum-frequency 
space in a form which is very appropriate for numerical calculations. 

We also perform numerical calculations on a lattice. Our numerical re- 
sults indicate that due to the decrease of the effective Coulomb repulsion 
at low frequencies when the quantum fluctuations are included, the state 
with an antiferromagnetic long-range order is more difficult to formed. We 
also compare the LO and NLO contributions to the fermion and boson self- 
energies. We find that the NLO contributions to the boson self-energy and 
the real part of the fermion self-energy are less than the LO ones, but as the 
Coulomb repulsion energy U is larger than about 2t, the NLO contributions 
can not be neglected. However, for the imaginary part of the fermion self- 
energy, it is found that the NLO contributions are comparable with the LO 
ones, even when the U is less than It. However, their signs are opposite and 
their sum almost approaches zero. 

Based on our calculations, we conclude that the 2PI effective action for- 
malism, which is popularly employed in the particle physics and field theories, 
is also an appropriate approach to describe the strongly correlated electron 
systems. Higher-order quantum fluctuations are easily included in this ap- 
proach. We must point out that one should be very careful about results 
obtained from the mean-field calculations of the strongly correlated electron 
systems, or even those including the LO quantum fluctuations. Because 
higher-order quantum fluctuations in the systems also play an important 
role, which can not be neglected. 
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